The electronic instabilities of the Iron-based superconductors: a variational 

Monte-Carlo study 
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We report the first variational Monte Carlo (VMC) study of the iron-based superconductors. 
We use realistic band structures, and the ordering instabilities/variational ansatzs are suggested 
by previous functional renormalization group calculations. We examine antiferromagnetism, super- 
conducting pairing, normal state fermi surface distortion and orbital order in the antiferromagnetic 
state. 
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The variational approach has a glorious history in con- 
densed matter physics. Examples of successful wave- 
functions include the BCS wave function, Laughlin wave 
function, Jastrow wave function, etc.. In the study 
of the cuprate superconductivity, Anderson's projected 
BCS wave function has been shown to capture many im- 
portant aspects of the cuprates[l|. Despite the successes, 
all variation approach have the draw back of being biased: 
i.e., it bases on the assumption that the few variational 
parameters built into the variational ansatz enable it to 
capture the essence of the true ground state wavefunc- 
tion. 

Understanding the pairing mechanism and possible 
electronic instabilities in iron-based superconductors 
have been a focus of interests in the past few years. This 
is not only for their high T c , but also for many similarities 
they share with the cuprates. Presently most of the works 
attribute the electron pairing to the antiferromagnetic 
fluctuations, and it is widely believed that the supercon- 
ducting (SC) gap function takes opposite signs on the 
electron and hole Fermi surfaces. This is partly because 
important consequences of the opposite pairing sign pre- 
dicted for the neutron@, |j| and STM exerimentsji, @ 
have now received experimental supportsJaLZf. 

On the theory side, except Ref.[8[ which assumes 
the iron-based compounds are doped Mott insula- 
tor, most of the studies are based on weak coupling 
approximations 19M20J . Experimentally there are reports 
favoring weak[21| and intermediate [22j electronic cor- 
relation in this class of compounds. Among Ref.[9|- 
|20j the functional renormalization group (FRG) [111] ap- 
proach has the virtue of being unbiased. It sums all one- 
loop particle-particle, particlc-holc diagrams as well as 
the vertex corrections. The results predict that in ad- 
dition to the AFM and SC instabilities, the iron-based 
compounds have propensity toward Fermi surface dis- 
tortion, magnetically coupled orbital order, and charge 
density wave orderings[18J. However given the fact that 



there is an evidence that the pnictides are intermediate 
coupling materials [22|, we must ask which, if any, of the 
weak coupling results in Ref.[9|420( are valid. One of the 
main purpose of this letter is to address the above ques- 
tion by performing a variational Monte-Carlo (VMC) cal- 
culation. In addition, while FRG only allows us to ac- 
cess the ordering tendencies, VMC allow us to quanti- 
tatively study the ordered state. We use the realistic 
bandstructurcs, and the variational ansatzs are guided 
by the predictions of Ref.[18j. We use partially projected 
wave functions where both the degree of projection and 
the order parameters are variational degrees of freedom. 
We shall focus on normal state Fermi surface distortion, 
AFM coupled orbital order, and superconducting pair- 
ing. (In FRG the tendency toward charge density wave 
order is the weakest.) 



The band structure we use is that of Kuroki et al. 15]. 
We model the electronic correlations by the Hubbard and 
Hunds types of local interactions. The Hamiltonian is 
given as follows: 
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where -ffhand and its parameters can be found in Ref. [151 ] 
and its Table I. In most of the paper we shall use 
(U 1 ,U 2 ,Jh) = (4, 2,0.7) eV. While these parameters 
arc compatible with Ref.[22j it is larger than what's re- 
ported in Rcf.[21|.The additional reasons for doing so 
are (1) similar large correlation parameters are con- 
cluded from LDA+DMFT calculations [23| although with 
controverserv|24J . (2) smaller interaction parameters will 
make the energy gain associated with the superconduct- 
ing and orbital order too small for our VMC accuracy. 
As the energy reduction due to normal state FS distor- 
tion is relatively larger, smaller interaction parameters, 
e.g. (Ui, U2, Jh) ~ (2, 1, 0.35) eV, were also used. 



We use the following partially-projected variational 
wave- functions, 



where 



IV-) = ^g^'Ptf. \tPmf) , (2) 

m v . (3) 






Pjv c fixes the total number of electron to N e , and < 
Si 1 52 < 1 suppress configurations according to the num- 
ber of electrons residing in the same and different on-site 
orbitals. The mean-field wavefunction \iPmf) depends on 
the type of electronic order we are studying, and it will 
be discussed in details later. 

The VMC was carried out on L 2 = 10 x 10 and 
18 x 18 lattices with periodic boundary condition. In 
choosing the doping level, we were careful to avoid de- 
generacy in the band energy between different Slater 
determinants. Standard Markovian chain Monte Carlo 
approach was used with Metropolis update algorithm. 
After an initial "equillibration" using 10 5 Monte-Carlo 
steps, about 10 5 ,10 6 and 10 7 samples were used in the 
study of Fermi surface distortion, orbital order and SC 
pairing respectively. The adjacent samples are separated 
by 5-L 2 Monte-Carlo steps with each showing an accep- 
tance ratio of about 0.23, which is enough to eliminate 
autocorrelation and thus guaranty an efficient sampling. 
Due to the higher demand of accuracy when studyin g SC 
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pairing, we have applied the "re-weighting" scheme 
The error-bars are estimated by calculating the variance 
of the energy expectation. 

Normal state FS distortion FRG predicts two lead- 
ing FS distortion: one preserves the 90° rotation symme- 
try and the other breaks it[lj|. The former shrinks both 
the electron and hole pockets and produces a relative en- 
ergy shift of the electron and the hole bands |26|. The 
later is the band version of orbital ordering, and is sug- 
gested to be stabilized by the AFM 18| ■ The mean-field 
state \iPmf) for the two types of FS distortion are the 
ground state of the following quadratic Hamiltonian 



#MF = X>k + x(k)W c 



(4) 
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where e£ is the bare band dispersion, n-^ aa 

and x (k) can be xo cos k x cos k y (90° preserving) or 

Xo (cos fez — cosky) (90° breaking). 

In the absence of AFM order, the results for 0.6% 
hole-doping on a 18 x 18 lattice using (Ui 7 U-2,Jh) = 
(4, 2,0.7)eV are shown in the main panel of Fig[TJ The 
black and red symbols represent E(xo) — E(0) as a func- 
tion of xo with gi_2 fixed at their optimized values. (The 
optimized (31,32) ~ (0.31(4), 0.76(4)) are nearly inde- 
pendent of Xi as the energy gain associated with their 
optimization is much larger than that associated with 
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FIG. 1: (Color online) Energy gain per site as functions of 
Xo for the normal state FS distortion, and the band and real 
space version of orbital ordering. The results were obtained 
on a 18 x 18 lattice under periodic boundary condition. The 
parameters gi, gi are fixed at their optimal values. The size of 
points represents the error-bar and the dotted lines are guide 
to the eyes. The main panel is for U\ — 4, U2 — 2, Jh — 0.7eV 
at 0.6% hole doping. In the inset the result for at 13% hole 
doping, with U\ — 2, U2 = 1, Jh = 0.35 eV and cos k x ■ cos k y 
form-factor is shown. 



optimization of X-) From Fig. [TJ it is clear that while 
both cos k x cos k y and cos k x — cos k y FS distortions gain 
energy, the former is the most energetically favorable one. 
This is consistent with the FRG prediction. The green 
symbols in FigfT] represent the energy reduction due to 
a real space version of orbital order (see later). We note 
that it is not energetically favorable in the non-magnetic 
state. At the optimal xo the total energy gain is about 
27meV per site. The negative xo moves the bands near 
r downward and the bands near M upward. As a result 
it shrinks both the electron and hole pockets (it actu- 
ally splits each electron pocket into two smaller ones). 
The distorted Fermi surface is compared with the undis- 
torted one in Fig. [2Ia,b). The above results are com- 
pletely consistent with the FRG prediction [18(. In the 
inset of Fig. [I] we show the E(xo) — E(0) versus xo plot 
for C/i = 2, C/ 2 = 1) Jh = 0.35 eV with cos k x -cos k y form- 
factor at 13% hole doping. The resulting relative energy 
shift of the bands near T and M is 90 — 126 meV which is 
in good quantitative agreement with the experimentally 
observed ~ 100 meV relative shift 26j, |29j. 

Orbital order in the AFM state The \iPmf) for the 
AFM ordered state is the ground state of the following 
mean-field Hamiltonian, 



Hmf — Ho 



Z^ \ GC i»o C ™o 



M^e^ R % (5) 



where H$ is given by either Eq. |4] or Eq. [6] 
Here Q = (tt, 0) is the AFM ordering wave vec- 
tor. As for M^v we used the mean-field result of 
Ref.[3fJ where the non-zero M^ are M^^,/! = 1, .., 5, 
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FIG. 2: The Fermi surfaces at 0.6% hole doping, (a) In the 
absence of FS distortion, (b) The distorted Fermi surfaces 
associated with x (k) = Xo cos k x cos k y . 
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FIG. 3: (Color online) The energy gain per site in the AFM 
state as function of xo for the three types of electronic orders 
studied in Fig. 1(a). The same lattice size and interaction 
parameters were used. The doping is 0.6%-electron doping. 
The variational parameters gi,(/2,M M „ are fixed at their op- 
timized values. The size of points represents the error-bar. 
The dotted lines are just guide to the eyes. 



and Mi5,M5i (= M15). The variational study is per- 
formed by keeping the ratio between the optimal mean- 
field parameters ,(M 22 , M 33 , M i4l M 55 , M 15 )/M n = 
(0.95, 1.03, 1.1, 0.76, -0.096), while allowing M n to vary. 
For (U u U 2 , Jh) = (4,2,0.7)eV at 0.6% electron- 
doping (on 18 x 18 lattice) we found (31,52) = 
(0.49(0), 0.76(6)), A/11 = lAeV, with a total energy re- 
duction of about 3.0eV per site, with nearly 2/x^ ordering 
moment. This value is significantly larger than the mea- 
sured ones for the stoichiometric compounds. The dis- 
crepancy can be due to the omission of the fluctuations 
in the orientation of magnetic moments and the ordering 
wavevectors , and/or the large values of the interaction 
parameters. 

In view of the strong atomic-like ordering moments in 
the AFM state obtained above, when studying the orbital 
ordering in the magnetic state we also adopt a real space 
version of orbital ordering[2j|, where \iPmf) is the ground 
state of 
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Here n iaa = c\ aa c iaa . We note that the orbital d xz ,d yz 
we use are 45° rotated from the dxz,dyz orbitals in 
Ref|15|. From Fig. 3, we conclude that in the presence of 
AFM, the real-space d xz /d yz orbital order is the most 
energetically favorable (we suspect this is due to the 
fact that the large, localized, ordering moment in the 
AFM state). It produces a total energy gain of about 
4.6mcV per site. At the optimal xo the occupation- 
number difference between the xz and yz orbitals is 
n<l xz — nd , = 0.20(1). This value is enhanced above the 
~ 0.15 occupation difference already present in the pure 
AFM state. Our result agrees qualitatively with a re- 
cent photoemission resul t [271 and a first principle Wan- 
nier function calculation. [28J . Comparing this result with 
the green symbols of FigJIJa), we conclude that this or- 
bital order is stabilized by the AFM. 

Superconducting pairing The \i/jmf} we use for 
studying pairing is the ground state of the following 
mean-field Hamiltonian 

Hmf = 5>£ - i»c)rw + Yl ( A k c k«t c -k Qi + h - c ) -(7) 

kacr ka 

We have studied four different types of gap functions: 



A£ = A (cos k. 
Ak = 



cos ky) , A£ = Arj cos k x ■ cos h 



A , A£ = A {cosk x — cos k y ) 
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In addition we only allow pairing in bands that cross 
the Fermi energy. This is an excellent approximation 
because, as shown below, the gaps are quite small. 

Our calculations for the SC state were carried out on 
a 10 x 10 lattice at 10% hole doping. The optimized 
{gii92) ~ (0.31(1), 0.76(2)). Again, to an excellent ap- 
proximation, the optimal value of 171.2 and [i c do not 
depends on A . The results for E(A ) - E(0) with 
{911 <?2, Mc) fixed at their optimal values are shown in 
Fig. |U The results suggest that while both conventional- 
s and d x 2_ y 2 pairing raise the energy, s± and extended-s 
lower it. Interestingly, our result suggests the extended-s 
(A (cosfc x + cos k y )) form factor is slightly favored (by 
0.1 ± 0.08 meV) over the s± (A cos k x ■ cos k y ) one. On 
the surface this contradicts the predictions of s± pairing 
form factor! However it is important to realize that the 
s± gap function is not synonymous to cos k x cos k y . In- 
deed, the s± form factor obtained from several weak cou- 
pling approaches ll|, [14J, [15| has strong variation around 
the electron Fermi surfaces. We believe the near degen- 
eracy of the extended-s and the cos k x cos k y form factors 
suggests the optimal pairing form factor is a linear com- 
bination of the two (hence is anisotropic on the electron 
pockets). However it is extremely computing time con- 
suming to verify this, and we have not been able to do 
it. In addition, we caution that the degree of the gap 
function variation on the electron pockets will depends 
on the values of the parameters. We cannot rule out 
that for other parameter sets the cos k x cos k y can be the 
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FIG. 4: (Color online) Energy gain per site as function of 
Ao for four different pairing form factors. The lattice size 
is 10 x 10, and 10% hole doping is studied. The other varia- 
tional parameters gi,g2, u c are fixed at their optimized values. 
Error-bars are shown. The dotted lines are guide to the eyes. 



leading pairing form factor. Finally, the optimal gap am- 
plitudes for the two symmetries are about 18meV and 
20 — 22meV respectively. Order of magnitude wise these 
values are not far from the experimentally measured gap 
scales. 

In conclusion, we have performed variational Monte- 
Carlo calculation to check the validity of our previous re- 
sults based on weak-coupling approximation. The results 
are qualitatively consistent. We caution that because of 
the time consuming nature of the calculation we are only 
able to study two sets of interaction parameters. Clearly 
we can not rule of the possibility that quantitative as- 
pects of the above results will be sensitive to the precise 
values of the interaction parameters. 
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